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Abstract 

The distribution of fragments produced in events involving the multifrag- 
mentation of excited sources is studied for peripheral Au + Au reactions at 
35 A.MeV. 

The Quasi-Projectile has been reconstructed from its de-excitation prod- 
ucts. An isotropic emission in its rest frame has been observed, indicating 
that an equilibrated system has been formed. The excitation energy of the 
Quasi-Projectile has been determined via calorimetry. 

A new event by event effective thermometer is proposed based on the 
energy balance. A peak in the energy fluctuations is observed related to the 
heat capacity suggesting that the system undergoes a liquid-gas type phase 
transition at an excitation energy ~ 5 A.MeV and a temperature 4 — 6 MeV, 
dependent on the freeze-out hypothesis. By analyzing different regions of the 
Campi-plot, the events associated with the liquid and gas phases as well as 
the critical region are thermodynamically characterized. 

The critical exponents, r, f3, 7, extracted from the high moments of the 
charge distribution are consistent with a liquid-gas type phase transition. 

PACS: 24.10.Pa, 25.70.Pq 



1 Introduction 

Nuclear multifragmentation and its possible connection to the occurrence of a phase tran- 
sition of the liquid-gas type has been the subject of intensive theoretical and experimental 
investigations j^, |2|, |3|, ^, |5|, |6|, |^. Theoretical studies indicate that infinite nuclear matter 
has an equation of state similar to that of a Van der Waals gas which is characterized 
by the existence of a liquid-gas phase transition. Moreover, recent experimental results 
show strong evidence for this occurrence in fragmenting nuclear systems. In particular 
the ^Z/74Z?/iV collaboration has measured a caloric curve resulting from the fragmentation 
of the Projectile and Target Spectators formed in peripheral collision Au + Au at 600 
and 1000 A.MeV, and found a behaviour expected for a first order liquid-gas phase tran- 
sition 1^, ^. The EOS collaboration has also measured a caloric curve and extracted the 
critical exponents of fragmenting nuclear systems produced in the collision of 1 A.GeV 
Au nuclei with a Carbon target, finding values consistent with the expected values at the 
critical point [pX|]. 

To experimentally study a nuclear system undergoing a phase transition events com- 
ing from the decay of a single source in a wide range of excitation energies have to be 
considered, so that it is possible to explore different kinds of de-excitation, from evapora- 
tion to multifragmentation towards vaporization. One can then hope to cross the critical 



region [g, 11 1 where the system manifests signals of a phase transition. 

In addition, the size of the studied source should be as large as possible to minimize 
finite size effects; and the collective radial flow (due to the compression of the nuclear 
matter in the first stages of the reaction) , as well as the rotational energy of the emitting 
system should be as small as possible. 

In central collisions, with reactions involving heavy projectiles and targets, it seems 
possible to overcome the difficulty concerning the size of the system. Indeed in Au + Au 
collisions at 35 A.MeV the Multics - Miniball collaboration has studied |12, [l^, |^, |l5| the 
characteristics of an equilibrated composite system with a size larger than 300 nucleons. 
Also in the Xe + Sn reactions the Indra collaboration |[l^, 17 1 examined an intermediate 



system with a size similar to the Au one. In these collisions, however, when the system 
is well identified, its excitation energy is distributed in a narrow range around the mean 
value so that it is impossible to build an excitation function with these data, except if 
the same collision is studied varying the beam energy in an almost continuous way. Even 



so, the collective energy contribution to the excitation energy depends strongly on the 
incident energy and it is not a trivial process to disentangle this component. 

Peripheral collisions seem to give a better chance to study the critical region, when 
selecting the Quasi-Projectile (QP) and Quasi- Target (QT) in the final state. In the reac- 
tion Au + Au at 35 A.MeV the maximum available centre of mass energy per nucleon is 
8.71 A.MeV which determines the scale of the QP and QT excitation energies. In principle 
one can observe sources with an excitation energy from a few hundred of A.keV to several 
units of A.MeV. Moreover, the radial collective energy is limited to ~ 1 A.MeV for the 



most central collisions |12], and should be smaller in peripheral and semi-peripheral colli- 
sions. The contribution of rotational energy can be evaluated in the multifragmentation 
region by comparisons with thermal predictions. 

Events associated to the Quasi-Projectile decay were already experimentally studied 



by looking for signals of a critical behaviour |18| , [1S| ] and anomalies in the caloric curve |20]. 
Several indications in favour of the occurrence of a critical behaviour were found by the 
study of the moments of the charge distributions [ p!8| , |l9|. They have been revealed mainly 
by the particular shape of the Campi scatter plot and by the presence of large fluctua- 



tions. Similar analyses were also performed |21] for peripheral collisions on the Xe + Sn 



50 A.MeV reaction. In the caloric curve analysis of Ref. [BOl, temperatures Tiso (calcu- 
lated from several double isotope ratios) increasing from 3.7 to 4.5 MeV with decreasing 
impact parameters were observed. These values resulted in agreement with Aladin and 
EOS data in the common range of excitation energy. 

Since the investigation of phase transitions implies the study of observables related to 



the critical behaviour on an event by event analysis p^ , ^ , ^, 25], the temperatures, as 
well as the excitation energies measured so far, cannot be used as the critical parameter 
as they are known only on the average. Indeed the experimental determination of these 
observables has been performed up to now in bins of other measured observables, such as 
Zbound [§, ^5 charged particle multiplicity p^, 26, 27, 28, 2£] and transverse energy [31|. 



In addition, the value of the measured temperatures could result from an average of 
different stages of the de-excitation [^, |2^. In the framework of the statistical model 
SMM [^, for instance, it has been shown that decaying systems with freeze-out tempera- 
tures of about 6 MeV give final isotopes with apparent double ratio temperatures on the 
order of 4 MeV, reasonably approaching the experimental values pO|. Only when consid- 



ering excitation energies much higher |^^ than those here considered the isotopic ratios 
give a good approximation of the freeze-out temperature. Other ambiguities arise from 
very recent experimental results on the temperature determination through the excited 
state populations |^. No effort is here made to investigate the reasons of the discrep- 
ancies between results obtained with different isotopic and excited states thermometers 
and to establish which are the most reliable ones. (For a discussion on this topic, see for 
instance Ref.s |^, 34, |35|]). 



In this paper we investigate in detail the disassembly of the QP formed in peripheral 
^^ Au + ^^^Au 35 A.MeV collisions. The data selection method guarantees in each event 
a unique emitting source of the fragments. 

A new effective event by event thermometer is proposed which allows for the first time 
the extraction of a "primary" caloric curve by backtracing the experimental information to 
the freeze-out configuration. It also allows for the computation of new thermodynamical 
observables, such as an effective heat capacity, through the energy fluctuations. A detailed 
analysis of the moments of the charge distribution is performed on an event by event basis. 



2 Experiment 



The experiment was performed at the National Superconducting Cyclotron Laboratory of 
the Michigan State University. Beams of Au ions at E/A = 35 A.MeV incident energy 
from the K1200 cyclotron were used to bombard Au foils of approximately 2.5-^ areal 
density. 

Light charged particles and fragments with charge up to the beam charge were detected 



at 9iab from 3° to 23° by the Multics array [36|, with an energy threshold of about 1.5 



A.MeV, nearly independent on fragment charge. Light charged particles, Z = \ and 
Z = 2 isotopes and fragments with charge up to Z = 20 were fully identified by 160 
phoswich detector elements of the MSU Miniball [37|, covering the angular range from 23° 



to 160°. The charge identification thresholds were about 2, 3, 4 A.MeV for Z = 3, 10, 18, 
respectively. 

The geometric acceptance of the combined array was greater than 87% of 47r. Dynam- 
ical calculations indicated [g^] that the apparatus is able to detect all the generated events 
up to a maximum impact parameter 6 ~ 11 fm, with the experimental triggering condi- 



tion (at least two fired detectors). Higher impact parameters 6 = 12 — 14 fm are detected 
with lower efficiency, being the Quasi- Target products below the energy thresholds and 
the Quasi-Projectile ones mainly flying at laboratory angles smaller than the minimum 
detection angle. 

We have studied in detail the efficiency for QP detection by performing calculations 
with the statistical SMM model [pO] for an Au source, excitation energies ranging up to 
9 A.MeV and freeze-out density po/3 (po stands for the normal nuclear density). For the 
velocity of the source in the laboratory frame, the value vqp was assumed from energy 
conservation: 



""QP = Y 2 {(^c.m. - f*Qp) + Vc.m. (1) 

where ec.m. is the available centre of mass energy per nucleon, Vc.m. the velocity of the centre 
of mass of the reaction in the laboratory reference frame and e^p is the QP excitation 
energy. 

The predictions were filtered through the software replica of the apparatus, taking into 
account the geometry, energy thresholds, Z-identification limitations and multiple hits. It 
resulted that for events where the total detected charge is at least 70% of the projectile 
charge, the event intrinsic efficiency for the detection of the QP products (ratio between 
the charged particle multiplicity detected by the apparatus and the one before the filter) is 
on the average 70% and that the ratio between the total charge bound in form of fragments 
after and before the filter is greater on the average than 93% (96%) for fragments with 
charge Z > 3 (6), respectively, over the whole range of the excitation energy e*. 

For these reasons the analysis presented in this paper is essentially based on the de- 
tected fragments. We will show in the next Section that, in our range of excitation energy, 
the limited efficiency for light-particle detection does not affect the calculation of observ- 
ables which characterize the decaying systems. 

3 Data Selection and QP characteristics 



As in the previous analyses of peripheral collisions [18, |19[ we selected the most peripheral 
collisions by requiring the velocity of the largest fragment in each event to be at least 75% 
of the beam velocity. After a shape analysis, i.e. an analysis of the variance of the fragment 
velocities |l^, ^ , the other IMFs of each event were considered as belonging to the QP 



if forward emitted in the ellipsoid reference frame. Only events where at least 70% of the 
Au charge was detected in the forward hemisphere of the ellipsoid were accepted for the 
QP source reconstruction. The velocity of the source was then calculated as the weighted 
mean of the IMF velocities. The last step is to add the light particles to the source. This 
has been done by doubling the contribution of light particles forward emitted in the source 
reference frame to minimize the contribution of pre-equilibrium particles p^, 00] . 

In Fig. 1 a), b) and c), we display the probability spectrum of cos{9fiow)i where Ojiow is 
the angle formed by the main eigenvector resulting from the shape analysis and the beam 
axis. Different panels correspond to different selections applied to the data. Fig. 1 a) 
corresponds to all the collected events, Fig.l b) to events with the constraint on the 
velocity of the largest fragment and finally Fig.l c) represents the cos{9fioyj) spectrum for 
the events satisfying both the conditions on the velocity of the largest fragment and on 
the total charge detected in the forward hemisphere. 

In Fig. 1 d), e) and f) a recently proposed observable |41, 42| is presented, which allows 
for the visualization of the source(s) in the analyzed events. The ensemble averaged charge 
density {pz{vpar)) is defined as: 

{Pz{v,ar)) = (^%^> (2) 

where pz{vpar) represents the event by event distribution in the velocity Vpar of the col- 
lected charge fraction. In Ref. E^] the preferred axis to project fragment velocity vectors 
was the main eigenvector of the ellipsoid reference frame, in our case the axis chosen is 
parallel to the QP velocity. The sums in the numerator and denominator of equation (|2[) 
have been limited to fragments {Z > 3). This is because light particles, symmetrized with 
respect to the QP velocity, do not contribute to the source identification. For our events 
this observable thus represents the distribution of the collected charge bound in fragments 
along the direction of the QP velocity. 

Panels d), e) and f) of Fig. 1 correspond to the same selections performed for the plots 
shown in Fig. la), b) and c), respectively. In all these panels a clear contribution from 
QP emission is seen, while the Quasi- Target products (below the energy thresholds) are 
lost, even if fast moving forward emitted light fragments can contribute to the observed 
distribution. The fragment emission in the mid-rapidity region decreases while applying 
the constraints on the largest fragment velocity and on the total charge detected in the 
forward hemisphere. By analyzing through gaussian fits |43| the pz distribution of Fig. 1 



f ) in different bins of impact parameter, we found indication that the possible pollution of 
fragments not belonging to the QP goes from about 2% for the most peripheral collisions 
to less than 10% for the most central ones. 

The distribution of the charge of the reconstructed QP (Zqp) is shown in Fig. 2 a). 
For the analysis presented in the following Sections we selected events where Zqp is 
within lb 10% of the Au charge (shaded area). For these events the mean measured total 
parallel momentum is 79% of the projectile linear momentum with a standard deviation 
7%. Fig. 2 b) represents the yield of different reduced impact parameters b/bmax after 
the selections leading to the plots of Fig. 1 c) and f) and the shaded area of Fig. 2 a). 
The ratio b/bmax has been calculated as in Ref. |Q through the detected charged particle 
multiplicity before the selection process except for a minimum bias trigger (at least two 
fired detectors). 

Being that the analyzed events are mainly peripheral with small contaminations from 
the mid-rapidity source and well detected from the point of view of the size of the QP 
and the linear momentum, the excitation energy has been calculated via calorimetry. 
Possible contributions of non thermal origin to the excitation energy have been evaluated 
by comparing the measured kinetic energies to thermal predictions. Indeed some pre- 
equilibrium light particles or fragments could have been included in the reconstructed 
Quasi-Projectile and contributed to its excitation energy. 

The excitation energy E* of the Quasi-Projectile in each event results from the energy 
balance between the initial stage at the freeze-out and the final stage of detected fragments: 

mo + E*= [J2imj + Ej) + ^(m„ + E^)] . (3) 

Here niQ is the mass of the QP, determined from the reconstructed charge assuming a 
charge-to-mass ratio as in the entrance channel, rrij, Ej {j = 1, . . . , Nc) are the masses and 
the kinetic energies (in the source reference frame) of the charged products (multiplicity 
-^c), rrin, En are the masses and the kinetic energies of neutrons. The rrij were calculated 
from the measured Zj through the numerical inversion of the Epax parameterization |p5[| , 
except for charges Z = 1, 2 where the measured values were used, mo, rrij, rUn take into 
account the mass excess. 

The number of free neutrons was obtained as the difference between the number of 
nucleons (Aq) of the QP and the sum of nucleons bound on the detected particles and 
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fragments ( A'^ = ^o — J2^j)- The average kinetic energy of neutrons was estimated with 
different approximations which lead to values of the excitation energy varying less than 
5%. These estimates were done: 



i) by assuming for £"„, [ p9| a Maxwell-Boltzmann thermal distribution, consistent with 
volume emission [^]. In equation @ X^-E-n = ^NnT = ^Nn\/E* /a has been as- 
sumed, where a = A/13 [E^ represents the level density parameter of the degenerate 
Fermi gas. 

ii) by assuming that in each event neutrons have the same mean kinetic energy as protons 
but decreased by the Coulomb barrier from the QP source. 

iii) by assuming, similarly to SMM model, that the Coulomb barrier depends on the ther- 
mal excitation energy of the emitting system (the light particle emission switches 
from surface evaporation to volume emission with increasing excitation energy). 
This requires applying a two-step calculation to the data. First we calculated the 
proton barrier with the estimate of E* described in i) and then, by solving equa- 
tion (H), we obtained the final value of E*. 

As already mentioned, possible contributions of non thermal origin were evaluated by a 
model comparison (see next Section). They resulted smaller than 10% in the whole range 
of the QP excitation energy. In the following, unless differently stated, we will use only 
the thermal part of the excitation energy i.e. the value obtained via calorimetry, reduced 
by the nonthermal contribution. We will show however, that in our excitation energy 
range, possible distortions introduced by the non thermal component do not change our 
main results and conclusions. 

Equation (^) was checked by applying it to SMM events and by comparing the obtained 
values to the input of the model. We found that calculated values differ from the input by 
no more than 3% both when applying eq. (P) to filtered or not filtered events. From this 
result we can conclude that, in our energy range, a high efficiency for fragment detection 
is essential for a correct determination of the excitation energy, while some inefficiency for 
light particle detection only slightly affects the calculation. 

In Fig. 2 c) we report the e* = E* /Aqp experimental distribution. The weight of the 
impact parameter on the experimental yield is evident. Therefore we will perform the 



study of the source characteristics and other observables in bins of e* (or h/hmax) narrow 
enough to safely neglect the weight of the impact parameter within each bin. 

In Fig. 2 d) the experimental correlation between the laboratory velocity of the QP and 
its excitation energy is shown. The line represents the values obtained from equation (|l|). 
The arrows represent the values of the beam velocity and of the centre of mass velocity 
in the laboratory reference frame. The experimental values are in agreement with the 
expected values, apart from some deviation at high excitation energy. Indeed our event 
selection could have partially eliminated events where the Quasi-Projectile was slowly 
moving with respect to the other fragment sources. 

To deeper investigate whether the QP can be considered as the unique source of the 
observed fragments, we have studied the charge density {pz) (equation (|2|)) in different bins 



of excitation energy (Fig. 3). We should to remark that recent theoretical studies |48] have 
shown that some deformation in the velocity spectra of light fragment can be ascribed to 
the Coulomb influence of the second source (QT) on the decay of the QP. This prediction 
could also explain the slight deformation of {pz) observed in Fig. 3 d). For the first three 
panels deviations of {pz) from a symmetric shape are negligible, since at low excitation 
energy the multiplicity of light fragments is low. At higher excitation energies (panel 
d)) the tail of the distribution is slightly stretched in the backward direction, as expected 
when light fragments minimize the potential energy through their position in the freeze-out 
volume. 

In conclusion the distributions shown in Fig. 1 and the QP characteristics shown in 
Fig. 2 and 3 make us confident that the utilized constraints do, in, select the events where 
the QP break-up is the dominant mechanism of fragment production. 

4 Comparison with a statistical model 

To verify the extent to which the reconstructed QP can be considered as the unique 
fragment source for the selected events, we compared the experimental Z-distributions 
with the predictions of a statistical model. 

The detailed description of the Statistical Multifragmentation Model (SMM) can be 
found in Ref . |g] . Here we only emphasize that the basic assumptions of the model are the 
equilibration of the decaying system and the statistical distribution of the probabilities of 
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the break-up channels. 

In Ref. [^ it was already shown that SMM can describe rather well the mean elemental 
event multiplicity N{Z) for Au + Au 35 A.MeV peripheral collisions. Here we want to 
look in more detail at the degree of reproduction of other observables. In principle, the 
detailed reproduction of the variances of the observables distributions should require the 



use of sophisticated backtracing procedures [15| which work in the multivariate space of 



the source characteristics. But this extraction of the sources distribution is beyond the 
aim of this paper. So we simplify our analysis by comparing experimental observables 
with SMM predictions coming from fixed values of the source size (^o = 197, Zq = 79, 
density /Oo/3) and for excitation energies continuously varied up to 9 A.MeV (8.71 A.MeV 
is the cm. available energy per nucleon for our reaction). For the velocity of the source 
in the laboratory frame the values from energy conservation were assumed (eq. ||) . 

In Fig. 4 we present the mean elemental event multiplicity distribution N{Z) for several 
bins of impact parameter. No relative normalization on the presented histograms was 
performed. The distributions have been only normalized to the total number of events 
both for data and SMM predictions. The solid line in Fig. 4 represents filtered SMM 
predictions, while the dashed one corresponds to not filtered ones. One can see that no 
distortions are introduced by the filter, since the apparatus has a quite high efficiency 
for QP products (well above the energy thresholds). With the considered QP ensemble 
the model reproduces quite well the experimental distributions, apart from the region of 
Fission Fragments at the most peripheral impact parameters. In particular it seems that 
the calculation gives less fissions and a more asymmetric mode with respect the data. 
However this is a particular channel existing at low excitation energy (about 2 A.MeV for 
the experimental data). Moreover the charge asymmetry between the two fission fragments 



can be infiuenced by a rotational motion of the emitting system [49| not considered in the 
present calculation. 

By removing both from data and model predictions the Fission Fragments (Fig. 5) 
it is clearly seen that the agreement of N{Z) is improved and that also the event charge 
partition, i.e. the distribution oi N[Z) for the heaviest fragments in each event, is very well 
reproduced. In Fig. 5 we show experimental events corresponding to the impact parameter 
range 0.8 < b/bmax < 0.9 but similar results have been obtained for all the b/bmax intervals 
shown in Fig. 4. 
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So we believe that this hmited lack of agreement between data and SMM is well 
understood and does not affect the multifragmentation region which is of interest for the 
analyses of this paper. 

As mentioned in the previous Section SMM predictions were also used, in each bin 
of impact parameter and for each charge, to evaluate possible contributions coming from 
pre-equilibrium light particles and fragments and/or from some collective energy. This was 
done by comparing the measured and predicted kinetic energies per nucleon. From Fig. 6 
it is evident that, while for heavy fragments the agreement between data and predictions 
is excellent, for fragments with charge smaller than 10 an extra kinetic energy is present 
in the data relative to the thermal+Coulomb contribution. These differences have been 
taken into account to evaluate the non thermal contributions to the calorimetric excitation 
energy. 

In Fig. 7 the experimental mean calorimetric excitation energy is compared with the 
SMM mean values (it the standard deviation) needed to reproduce the experimental charge 
distribution and the charge partition. For this comparison we eliminated fission events 
both from data and predictions. For all the bins of the impact parameter a very good 
agreement is evident either including or excluding the collective energy. 

Our experimental partitions show a high degree of thermalization, apart from a small 
contribution of extra-kinetic energy for light charged particles and small fragments. Suc- 
cessive applications of thermalization hypothesis to the data should take into account that 
some nucleons and light particles can have escaped the system through surface emission 
before the thermalization. Even if this effect should be more important at excitation 
energies higher than those here considered in our following analysis it will be taken into 
account. 

The high degree of reproduction of the observables in each bin of impact parameter 
shows that SMM gives an accurate description of the phase space experimentally ex- 
plored [12, 15, 17, |5^, 51]. These findings make the analysis of temperatures even more 



intriguing. 
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5 Temperature of the QP source 

As we have shown the relevance of the statistical approach to describe the selected events, 
we will now present an analysis about thermodynamical observables of the multifragment- 
ing source. In particular we will extract from the data information about the freeze-out 
temperature of the decaying system. This experimental temperature will be compared 
with the predictions of two statistical models MMMC 0] and SMM [^, which use some- 
what different assumptions on the freeze-out configurations. 

The first model [0] assumes that the break-up volume is fixed for all partitions. Pri- 
mary fragments are produced in ground and low excited states and they de-excite only 
by evaporating neutrons. Their de-excitation is treated in a microcanonical way inside 
the freeze-out volume. In the second model [^ one assumes that the break-up volume 
grows with the number of the fragments. Primary fragments have the same temperature 
as the whole system and in a later stage, during the Coulomb acceleration, they de-excite 
by evaporating neutrons and charged particles or by secondary Fermi break-up. In ad- 
dition, at low excitation energies (typically up to 4-5 A.MeV), some final fragments are 
sequentially emitted in the decay of the heavy residual nucleus present at freeze-out. 

From a general point of view, one should bear in mind that partitions at freeze-out are 
different from the asymptotic measured ones. Therefore, following the two different freeze- 
out assumptions of the models we reconstructed primary partitions and we estimated the 
temperature Q of the QP in each event looking for a new thermometer mainly based on 
fragments. The energy balance at the freeze-out time reads: 



mo + E* 



^ m, + ^ Er' + E,oui + ^(M - l)e 



(4) 



2 

Here tjiq is the mass of QP, E* is the excitation energy of the source calculated via 
calorimetry through equation (y), m, (i = 1, . . . , M) are the masses of primary fragments, 
light particles and neutrons, M is the total multiplicity, i?*"* = a@^ is the internal energy 
of primary fragments (Z > 3), E^oui is the Coulomb energy of the partition calculated by 
randomly positioning non overlapping primary fragments in the freeze-out volume. The 
free parameters entering the equation at the freeze-out are the level density a and the 
freeze-out density p. 

To evaluate for the data the different terms of eq. (Q) we followed the two different 
assumptions of the models MMMC and SMM Q. To mimic the MMMC assumption 
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we assumed for primary fragments (Z > 3) a charge-to-mass ratio as in the entrance 
channel. In the SMM hypothesis we shared among final fragments the total charge detected 
in form of light particles and we assumed, as before, a charge-to-mass ratio as in the 
entrance channel. This sharing was performed in two different approximations. In the 
first it was assumed that primary fragments will evaporate particles proportionally to 
their charge, while in the second approximation the same amount of charge (detected in 
form of light particles) was distributed to all the fragments to mimic a higher Fermi-break- 
up probability. These two assumptions give results in agreement within few percent. In 
addition, following Ref. [y] we allowed compound nucleus emission to simulate a secondary 
emission for the detected fragments. We have bound the lightest fragment in each event 
with the largest one with a probability 



■W/<^? 



P(e*) = e-2i^ /-r. (5) 

The width of the gaussian P{e*) was assumed to be cr = 4 A.MeV, based on the behaviour 
of the experimental fragment multiplicity as a function of the excitation energy. Indeed 
the two-fragment channel is experimentally seen in a wide range of excitation energies, 
up to ~ 4 — 5 A.MeV, until the largest fragment in each event becomes smaller than the 
Fission Fragments. 

defined algorithmically by equation (^ has to be interpreted as an event by event 
estimator of the microcanonical temperature M, M. For each E*, Q has a spread accord- 
ing to the kinematical properties of the partitions. Similarly to SMM model, indeed, the 
partition temperature depends on the amount of kinetic energy still available after the sub- 
traction of the Q-value to create a partition and of the Coulomb energy of the interacting 
fragments. Averaging Q over all the partitions at fixed E* one gets an estimate of the 
microcanonical temperature of the system. 

Equation (|^ was first applied to the model events to estimate how this backtracing 
from t = oo (final partitions) to t = (freeze-out configuration) is able to approximate the 
original correlation between the microcanonical temperature and the excitation energy. 
The method resulted to work quite well as it reproduces the microcanonical tempera- 
ture within 10% in the case of MMMC model and within 5% in the SMM case. Since 
equation (^ resulted reliable, hereafter we use the notation T instead than (0). 

In Fig. 8 the experimental caloric curve, calculated through the total energy balance 
(eq. @), is shown for the two freeze-out hypotheses SMM-like and MMMC-like. A freeze- 
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out density p = po/3 has been assumed for the first hypothesis on primary fragments 
and a density p = po/G for the second one. Concerning the level density parameter we 
used a = f{A), with / continuously varying from A/ 8 for the lightest nucleus A = 7 
to A/12 for the heaviest one A = 197 to take into account different surface and volume 
contributions Q. In this figure we also report the caloric curves of SMM and MMMC 
models, calculated for filtered events. The effect of the experimental acceptance on the 
caloric curve (an increase of the order of 15% at the highest excitation energy) is to 
suppress events with a large (> 7) number of fragments which correspond in the model to 
partition temperatures lower than those corresponding to a limited number of fragments. 
Therefore the expected plateau-like behaviour |0, ^ is no more evident. 

Despite the difference in the absolute values of the temperature many similarities exist 
between the experimental (T, e*) reconstructed with the two hypotheses SMM-like and 
MMMC-like. In both cases at small excitation energies T oc ^/e* and becomes flatter 
for increasing energies even if the energy of this transition depends on the freeze-out 
hypothesis. If we consider the behaviour T oc \/e* as an indicator of a liquid phase, we 
can interpret this change as a first qualitative sign of a phase transition. 

To check the influence of the collective energy present in the experimental kinetic 
energies and of the parameters a and p used in equation (0), we have compared the ex- 
perimental caloric curves obtained under different conditions. We have found that any 
combination of parameters does not modify significantly the experimental (T, e*) corre- 
lations, the main change on the value of the freeze-out temperature arising from the 
freeze-out reconstruction hypothesis. 

It is interesting to note that the hypothesis SMM-like applied to the data gives values 
of the temperature systematically higher than those given by the hypothesis MMMC-like. 
This can be explained by the fact that, when applying eq. (Q), we always or never bind the 
light charged particles to final fragments to get primary fragments. Since however it is not 
experimentally established which is the excitation energy where a transition between the 
two regimes could take place we can consider our experimental caloric curves as an upper 
and a lower limit to bracket the freeze-out temperature and we estimate the accuracy of 
our reconstruction to be about 2 MeV. The contribution of secondary decay to the light 
charged particles and fragments yield has therefore to be evaluated before a more precise 
determination of the freeze-out temperature can be done. 
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The fact that SMM and MMMC for a given excitation energy and density produce 
approximately the same asymptotic partitions |^, ^ suggests that static variables do not 
allow to discriminate between the two freeze-out hypotheses. Therefore we investigated 
dynamical observables that could keep a memory of the different fragment production 
mechanisms, such as for instance the correlation functions of the reduced velocity: Vred = 

\Vi-Vj\ 



It is well established [13, 54 1 that at small values of the relative velocity the two- 
fragment correlation functions of the reduced velocity are quite sensitive to the interfrag- 
ment emission time and to the density of the emitting system. Here we want to exploit 
another property of this kinematical region, called Coulomb hole. Namely a partial filling 
of the Coulomb hole would reveal a decrease of correlation among fragments due to their 
late decay through light charged particle emission. 

We show in Fig. 9 the correlation functions of the reduced velocity for fragments with 
charge 4 < Z < 15 for theoretically calculated events at excitation energies 4 A.MeV and 
6 A.MeV. In the upper panel MMMC predictions are presented and in the lower panel 
SMM events are analyzed. In this case the codes were run at the same density (po/6) to 
enhance in the Coulomb hole region only effects due to secondary evaporation. 

The sensitivity of this dynamical observable to the different hypotheses of fragment 
production is evident. We clearly see how the SMM assumptions about Fermi break-up 
and light particle evaporation decreases the correlation among fragments with respect to 
the MMMC assumption of fragments which freeze only through neutron emission. 

Since the statistics of the measured events presented in this paper is too low to perform 
a calculation of correlation functions in different bins of the excitation energy, we present in 
Fig. 10 only two large intervals of e* which should be representative of a possible change 
in the fragment production mechanism. In the upper panel of Fig. 10 the correlation 
functions for e* < 3.5 A.MeV are shown and in the lower panel the same observable for 
5 < e* < 7 A.MeV. Both panels refer to fragments with charge 4 < Z < 15. 

Even if the statistics is poor, the difference between the two correlation functions is 
evident. The width of the Coulomb hole increases with increasing e* indicating a transition 
to a prompter emission of fragments. 

The experimental correlation functions of the reduced velocity cannot be directly com- 
pared with the predicted ones presented in Fig. 9. Indeed the model calculations are done 
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at a fixed excitation energy, while in the data the observables result from an average of 
events with different excitation energies. In addition an investigation about the density of 
the decaying system should be performed before comparing dynamic observables ||T^, but 
this is beyond the aim of this paper. We can only claim that data show a lower correlation 
among fragments than the one seen in MMMC model. The hypothesis of hot fragments 
with an important secondary emission seems more reliable, at least for our dataset. We 
should note that in Ref. p^] a sophisticated correlation analysis suggests the existence, in 
the multifragmentation of highly excited sources, of hot primary fragments emitting light 
charged particles. 

More statistics is needed to perform a quantitative comparison between experimental 
and predicted correlation functions in bins of excitation energy. We hope that measure- 
ments with high statistics will be soon performed by new Air detectors ^ with high mass and 
energy resolution, with the aim of studying, for a well defined source, the transition from 
sequential to prompt fragment emission through the correlation functions of the reduced 
velocity. 

To summarize this Section, we have shown that different freeze-out assumptions lead to 
different back-traced freeze-out temperatures and thus to a different origin of the measured 
partitions and of the isotopic yields. From correlation functions we got the indication that 
dynamic observables could be better reproduced by models where "hot" primary frag- 
ments are generated at the freeze-out, decaying later through light particle and fragment 
emission. In our range of excitation energies the isotopic temperature seems to be a mea- 
sure of the mean temperature of the emitting systems at a stage of the reaction later than 
the freeze-out. 

6 Signals of phase transition 
6.1 Fluctuations 

After the study of average temperatures, we analyzed fluctuations of (0, e*). It is hopeful 
indeed that Q fluctuations from event to event will provide a deeper understanding of the 
thermodynamical properties of the decaying systems. In Fig. 11 we show the (0, e*) scatter 



■'For instance the detector Chimera |Q, which will start to be operating in 1999 at the Super- 
conducting Cyclotron of the Laboratorio Nazionale del Sud, Catania 
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plots for data and SMM events, calculated in a larger energy range. It is clearly seen that 
in both cases the distribution has a broadening in the region around 4-^6 A.MeV. 

One of the observables related to energy fluctuations is the specific heat capacity. A 
peak for this observable at the phase transition was anticipated by theoreticians more than 
10 years ago [^, ^ and more recently studied in the framework of different models |58, 



1,0 



For a canonical ensemble the specific heat capacity at constant volume can be calcu- 
lated as: 

^v=j^ = ^ • (6) 

Since it is expected that Cy shows a peak at the phase transition temperature, this gives a 
hint that also in our case a'^/T'^ may be relevant to provide information in addition to the 
average characteristics. To better investigate this point, we study the energy fluctuations 
extracted from data on an event by event basis. 

We should nevertheless emphasize that in order to study energy variances, in our 
microcanonical type of framework, we have to consider fixed bins of Q and then combine 
various partitions with different energies. In this way we deal with energy fluctuations 
even if these fluctuations do not strictly coincide with those of a canonical ensemble. The 
absolute value of the total energy variance a'^/T'^ cannot in our case be directly interpreted 
as a heat capacity but it will rather be related in a non linear way both to the heat capacity 
and the latent heat of the transition. In any case even for a microcanonical system the 
location of the peak, if any, will sign the occurrence of a phase transition |]60t| . 

In Fig. 12 we plot the experimental a'^/T'^ as a function of the temperature for the 
two freeze-out hypotheses and for the different values of the parameter a used in the 
calculations. We show also calculations where the extra-kinetic energy was subtracted 
or included in the fragment kinetic energies. The dashed lines in Fig. 12 correspond to 
regions where the statistical error of the events in the bin of temperature is larger than 
10%. In all the cases al^/T"^ shows a peak and this qualitative behaviour is not modified 
if we vary the parameters /o, a or we include the extra-kinetic energy. 

One may worry if the decrease of fluctuations at high temperature is not trivially 
induced by the shape of the experimental total energy distribution: a monotonically in- 
creasing function could be biased by the decreasing weight of high excitation energies. 
To explore this possibility we have artificially flattened the excitation energy distribution 



giving to the detected events an energy dependent weight for the calculation of eq. (|6|). 
The result is a deformation of the shape of the fluctuations without any shift in the loca- 
tion of the peak. To complete this analysis and to verify if experimental inefficiencies and 
limitation at high excitation energy could modify the distribution of this observable, we 
report in Fig. 13 o"|;/T^ calculated for SMM not filtered events under different conditions. 
The lack of events at high excitation energies (temperatures) modifies the upper tail of 
the distribution and reduces the height of the maximum but again does not produce any 
shift in the position of the peak which appears to be robust. 

We can then conclude that the observation of a peak in the experimental effective 
specific heat, defined by eq. (g) and shown in Fig. 12, reveals a phase transition for the 
studied system. The two different freeze-out hypotheses lead to transition temperatures 
differing by about 1.5 MeV but almost identical excitation energies. However, one should 
remember that the ensemble of experimental events is limited in excitation energy. Data 
at a higher beam energy, leading to a wider excitation energy distribution, would be most 
welcome to confirm this finding. 

Another way to characterize fluctuations is to analyze the behaviour of (0, e*) for dif- 
ferent fragment partitions. One of the most powerful methods to characterize the critical 
behaviour of a system undergoing multifragmentation is the method of conditional mo- 
ments introduced by Campi pH]. The moments of asymptotic cluster charge distributions 



are defined as: 



and the normalized moments as: 



mu 



Y^Z^n^^\z) (7) 

z 



4^^=<V-1^^ (8) 

where n^^'{Z) is the multiplicity of clusters of charge Z in the event j, and the summation 
is over all the fragments in the event except the heaviest one which corresponds to the bulk 
liquid in an infinite system. 

Fig. 14 a) shows, for each event j, the experimental scatter plot of the logarithm 
of the charge of the largest fragment {ln{Z^^)) as a function of the logarithm of the 
corresponding second moment {ln{m2 )) (Campi scatter plot). If the system experiences 
a phase transition this plot should exhibit two branches: an upper branch with an average 
negative slope corresponding to under-critical events and a lower branch with a positive 
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slope that corresponds to over-critical events. The two branches should meet close to the 
critical point of the phase transition [22, 24, 25 1. 



As already done in previous studies of peripheral collisions [18, 19 [, to characterize the 
two branches and the intermediate zone, we have made three cuts in this plot selecting 
the upper branch (Cut 1), the lower branch (Cut 3) and the central region (Cut 2) and 
analysed the events falling in each of the three zones. It was already shown [l^, ll^ that 
the fragment charge distributions obtained in the three zones exhibit shapes going from 
a "U" shape in Cut 1, characteristic of evaporation events, to an exponential one for 
Cut 3, characteristic of highly excited systems. In Cut 2 a power law fragment charge 
distribution Z~'^ (r ps 2.2) was observed as expected according to the Fisher's droplet 



model for fragment formation near the critical point of a liquid-gas phase transition [61[. 

Here we argue that the Campi plot can be used also for a calorimetric analysis, to study 
the thermodynamical properties of the emitting Quasi-Projectiles for events belonging to 
the different cuts even if each selected region is representative but not exhaustive of a 
particular class of events. 

The upper branch contains events having an average reduced impact parameter 0.98, 
mean excitation energy 1.4 A.MeV, while the lower branch consists of events having 
{b/bmax) = 0.67, (e*) = 6.4 A.MeV. The region where the two branches meet is populated 
by events with (b/bmax) = 0.81, (e*) = 4.3 A.MeV. Remark that b/b^ax ~ 0.8 was found 
in theoretical calculations based on Classical Molecular Dynamics [p3| as the value of the 
impact parameter where the QP shows the signals of a critical behaviour. The temperature 
distributions of the QP depend on the assumptions made for the freeze-out configuration 
as discussed in the previous Section. Following the hypothesis where both neutrons and 
light charged particles are emitted by primary fragments we obtain for the cuts 1, 2, 3 
on the Campi scatter plot mean partition temperatures of 3.9, 5.8, 6.4 MeV, respectively. 
Following the hypothesis of fragment de-excitation only through neutron emission we 
obtain for the three cuts mean partition temperatures 2.9, 4.0, 4.2 MeV, respectively. 

The contributions to the caloric curve in the three cuts are shown by diff'erent colours in 
Fig. 14 b) for experimental events. For comparison in Fig 14 c), d) we show the Campi plot 
and the caloric curve for SMM filtered events. It is clearly seen, both for data and SMM 
events, that for cuts 1 and 3 there are unusual events (although with low probability) which 
lie far from the average T(e*) behaviour. These are unusual compound-like states with very 
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large temperatures and multi-fragments events with low temperatures. For these events 
one could do the analogy respectively with an overheated liquid and a super-cooled gas 
for ordinary liquid-gas phase transition. Their yield in the present experiment is too small 
to draw reliable quantitative conclusions. Further studies on B fluctuations will provide a 
deeper understanding |6C] of the thermodynamical properties of the under-critical, critical 



and over-critical events of the phase transition in finite systems. 

6.2 Critical exponents 

If multifragmentation is a manifestation of a phase transition in nuclear matter one can 
try to extract from the data a set of critical exponents characterizing this phase transition. 
Strictly speaking critical exponents are defined for 2-nd order phase transitions where the 
correlation length diverges. Indeed a liquid-gas phase transition is 2-nd order only at 
the critical point. Nevertheless, as we have seen in Fig. 11, a large transition region with 
wide energy fluctuations is populated by the selected multifragmentation events. One may 
hope that some events fall in the vicinity of the critical point. On the other hand in finite 
systems the universal critical behaviour can be also induced by finite size effects (see e.g. 
lattice based calculations |^]). 

For systems undergoing a second order phase transition the moments of the charge 
distribution (eq.s (^) and (0)) mi. and S^ should reveal some specific correlations p^ and 
in particular for k > 2 they should diverge at the critical point or show a peak, due to 



finite size effects. In previous papers |18, |l^, 32| peaks were found for the second moment 



1712, the reduced variance 72 and the normalized variance of the largest fragment ajsry as a 
function of the charge particle multiplicity. These are indications that a critical behaviour 
could occur. 

Critical exponents |1C, ^ are characteristic of the specific universality class of the 



transition. Here we recall that the exponent r gives the slope of the power-law distribution 
at the critical point, the exponent (3 fixes the behaviour of the largest cluster close to the 
critical point and finally the exponent 7 governs the strength of the singularity of m2 at 
this point. These quantities are connected by universal relations [^, |6^ . 

We could not extract 7 from the second moment of the charge distribution, as per- 



formed in Ref.s |1C, |62], because in our range of excitation energies we have too few 



events in the over-critical region to constrain the value of 7. Instead of this the ratio 
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/3/7 has been calculated by fitting the correlation ln{Zi,ig) versus ln{{S2)), mainly in the 
"liquid" region. The value of r has been obtained from the correlation between moments 
with different order S^ and Sj (l > k,k > 2), since the slope of the log-log correlation is 
(r — / — l)/(r — A; — 1). The exponent r has been also calculated by fitting the Z-distribution 
in the "critical" region (Cut 2 in the Campi-scatter plot). The exponent (3 has been ob- 
tained by the fit of the mean value of the charge of the heaviest fragment Z{,ig as a function 
of the excitation energy (temperature), assumed to be the critical parameter. 

Before calculating (3/^ from ln{Z},ig) versus ln{{S2)), great care was payed to the 



elimination of fission events, following the warnings given in Ref. [24|. These events indeed, 
enhancing the value of the moments, do not allow a proper determination of the critical 
exponents. From this study it resulted that there are two regions not affected by the 
criterion chosen to eliminate fission events: the higher part of the liquid branch and the 
lower branch. These regions have been considered to determine /3/7. By fitting (Fig. 15 a) 
the upper part of the correlation we obtain (3/'j = 0.29 it 0.01, while from the lower branch 
we get 1 + /3/7 = 1.31 it 0.07. In the latter case the error is large, because the fitted region 
in our case is limited (we do not observe vaporization events). 

We have calculated the ratio f3/^ also for SMM predictions. For not filtered events we 
get /5/7 = 0.28 ib 0.01 for the upper branch while from the lower branch, where at high 
excitation energy the model gives mainly vapor, we obtain 1 + P/j = 1.28 it 0.02. For 
filtered events we obtain f3/-f = 0.28 ± 0.02 and 1 + (3/j = 1.28 ± 0.03. It is important 
to note that the fits of the two branches give the same values of P/'j even in the case of 
filtered events and that these values are in agreement with those found from data. The 
theoretical value expected for a liquid-gas phase transition is P/'j = 0.265 it 0.006. 

The exponent r has been calculated from the correlations S3 - S2 and 5*4 - ^3. In 
the case of SMM events the fit of 5*3 — 5*2 gives r = 2.168 it 0.002. Experimentally, by 
fitting S3 versus S2 (Fig. 15 b)) we obtain r = 2.12 it 0.02. To verify the extent to which 
the experimental inefficiencies for light particle detection influence the value of r we have 
calculated the same correlation by correcting for the mean efficiency of the apparatus 
to the contribution of the light particles. With this correction we get r = 2.13 it 0.02 
compatible with the previous value. By fitting higher moments, like S4 versus S3, we 
obtain t = 2.13 ±0.01. 

In Fig. 15 c) we report the correlation 5*3 — S2 for the events falling in the regions 1, 
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2 and 3 of the Campi scatter plot (respectively squares, points, open circles). We observe 
that the events falling in the region 2 cover the whole range of the 5*3 — S2 correlation 
thus determining r. Evaporation events belonging to the Cut 1 fall in the region of small 
values of 5*3 and ^2 and multifragmentation events belonging to the Cut 3 fall in the upper 
corner of the correlation. 

We show in Fig. 15 d) the charge distribution for events belonging to the Cut 2 of 
the Campi-scatter plot. The line represents the power-law fit of the distribution which 
gives T = 2.13 lb 0.08, in agreement with the value 2.22 it 0.01 obtained for r from the 



relation r = 2 + 1/(1 + 7//3). It was claimed by several theoretical works |11, 64] that it 
is difficult to compare the exponents obtained for finite systems with the expected ones 
corresponding to the relevant universality class, due to the deformations induced by finite 
size effects. However the slope P/j [^, as well as the related exponent r [^, rapidly 
approach the value expected in the thermodynamic limit for increasing size of the system. 
As previously mentioned, to extract from data another critical exponent, we need to 



make an assumption about the critical parameter. It is predicted [^, 25] that Z^ig, the 
bulk in the liquid region, behaves as: 

Zbig ~ e" (9) 

where e (> 0) is a variable characterizing the vicinity to the critical point. In percolation 
€ = p — pc (p > Pc) is the deviation from the threshold probability Po while in thermal 
transitions of the liquid-gas type in the canonical ensemble e = T^ — T [T < T^) is the 
deviation from the critical temperature. 

Since experimentally e is not directly observable, a monotonic function of e can allow to 



determine critical exponents. For instance the multiplicity of produced cluster j|l^, 25, 62] 
has been often used. An alternative possibility is the thermal excitation energy e* which 
could be the proper thermal critical parameter in the micro canonical ensemble ]p4| , |65[] . 
This observable, together with the estimated temperature extracted from equation (|^, 
has been used in this work for the determination of the exponent /?. 

In Fig. 16 a) the correlation Z^jg — e* is shown, together with the fit resulting from 
equation (|). From this fit we obtain e*^^^ = 4.5 ± 0.2 A.MeV and {3 = 0.33 ± 0.04. We 
want to emphasize that in the range of excitation energy used for the fit (1 — 4 A.MeV), 
the possible contribution of collective energy to the thermal excitation one is negligible so 
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that the value of e*^jj is actually reliable. In Fig. 16 b) we show another representation of 
the behaviour of the size of the largest fragment as a function of the critical parameter: 
In^Zbig) vs. ^?i-(e*rjt — £*)• This representation should enhance the deviations from the 
power law behaviour (eq.( ^)). 

If instead of e* we use \/e* as critical parameter (which roughly corresponds to a 
temperature in the liquid side) we obtain e*^jj = 5.0 it 0.3 A.MeV and /? = 0.33 it 0.02. 
By using the reconstructed temperature as critical parameter (Fig. 16 c) and d)) we get 
Tcrit = 6.0 ± OAMeV, (3 = 0.35 ± 0.07 with the freeze-out hypothesis SMM-hke and 
Tcrit = 4.6 ± 0.3 MeV, (3 = 0.30 ± 0.07 in the MMMC-like hypothesis. 

The value of (5 extracted from SMM predictions is /5 = 0.344 it 0.004, in agreement 
with experimental data. 

The fact that a power law behaviour with the same (3 exponent is found for Z^yig when 
working with e* or ^/e* or T should not be surprising since this behaviour is expected 
only in the close proximity of the transition point. However this result confirms the fully 
statistical nature of the data, the coherence of our analysis and the reliability of eq. (Q) 
to calculate the freeze-out temperature. 

In conclusion we have shown that the correlation between the largest cluster and the 
excitation energy or the temperature leads to values of the critical parameter fulfilling our 
hypothesis on the relationship between e* and T. The values of e*^jj and Tent correspond 
to the mean values obtained in the region 2 of the Campi-scatter plot. These values 
correspond to the transition region estimated from the analysis of the fluctuations, where 
the effective specific heat shows a peak. The values of /3/7, j3 and r are compatible with 
theoretical predictions for systems undergoing a phase transition. 

We want however to remark that the extraction of critical exponents has been mainly 
based on experimental correlations related to the liquid branch of the Campi scatter plot. 
We think that an analysis similar to the one presented in this paper, performed also for the 
over-critical branch of the Campi scatter plot, could give a more complete understanding 
of the nuclear phase transition. 
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7 Summary and Conclusions 

In this paper we have presented a comprehensive study on the Quasi-Projectile multifrag- 
mentation in peripheral collisions for the reaction Au + Au at 35 A.MeV incident energy. 
The data selection constraints and the source reconstruction technique have been chosen 
to minimize preequilibrium and/or neck contamination in the predominantly binary events 
detected. 

The high degree of equilibration attained by the decaying QP system is also demon- 
strated by the agreement with statistical multifragmentation model predictions. The exci- 
tation energy and temperature of the source have been evaluated through a detailed energy 
balance at the freeze-out time. On an event by event basis the correlation between the 
temperature and the excitation energy of the system has been calculated by backtracing 
the experimental information to the freeze-out time. 

From an experimental point of view it has been shown the relevant importance of the 
fluctuations to characterize the multifragmentation phenomenon. A peak of the "canoni- 
cal" heat capacity indicates that the system undergoes a liquid-gas type phase transition 
at an excitation energy ~ 5 A.MeV and a critical temperature 4 — 6 MeV, dependent on 
the freeze-out hypothesis. 

Critical exponents, extracted from the moments of the charge distribution in the liquid 
part of the experimental Campi plot, have values compatible with theoretical predictions 
for systems undergoing a phase transition. 
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Fig. 1: 

Probability distribution of cos{9fiauj) (upper panels) and distribution of charge density 
{pz{vpar)) (lower panels) as a function of the fragment velocity along an axis parallel to 
the QP velocity. The arrow indicates the mid-rapidity velocity. 

Panels d), e) and f) correspond to the same selections performed for the panels a), b) and 
c), respectively (see text). 
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Fig. 2: 

a) Charge distribution of the QP. The shaded area (70 < Zqp < 88) represents the events 
analyzed in the present work. 

b) Yield of the reduced impact parameter. 

c) Experimental distribution of the QP excitation energy, obtained with eq. (^). 

d) Experimental correlation between the laboratory velocity of the QP and its excitation 
energy. The dashed line represents the values obtained from eq. (|l|). 

Panels b), c) and d) refer to events belonging to the the shaded area of panel a). 
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Fig. 3: 

Charge density {pz{vpar)) distribution as a function of the fragment velocity, along the 
axis parallel to the QP velocity. Panels a), b), c) and d) refer to e* = 1 — 2 A.MeV, 
e* = 2 - 3 A.MeV, e* = 3 - 6 A.MeV and e* =6-8 A.MeV, respectively The 
mid-rapidity source is centered at ~ —4 cm/ns. 
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Fig. 4: 

Mean elemental event multiplicity N{Z) in different b/bmax intervals. The circles show 
experimental data and the solid/dashed histograms the results of SMM filtered/not filtered 
predictions. 
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Fig. 5: 

Mean elemental event multiplicity N{Z) and distribution of the three heaviest fragments 
in each event, ordered for decreasing size (Zi > Z2 > Z3) for 0.8 < b/bmax < 0.9 with 
fission events removed. The circles show experimental data and the lines the filtered SMM 
predictions. 
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Fig. 6: 

Measured (circles) and predicted (lines) mean kinetic energy per nucleon as a function of 
the fragment charge for different b/bmax intervals. 
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Fig. 7: 

Mean excitation energy per nucleon and its standard deviation as a function of the reduced 
impact parameter. Points show the experimental mean values (+ the standard deviation), 
squares are the excitation energy decreased by the collective component (- the standard 
deviation). Solid/dotted lines are filtered/not filtered SMM mean values it the standard 
deviation. Fission events were eliminated from both data and predictions. 
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Fig. 8: 

Experimental temperature-excitation energy correlation. Circles represent the average 
values calculated with equations (y) and (^ in the SMM-like freeze-out hypothesis; squares 
in the {MMMC-like) hypothesis. Dashed and dot-dashed lines represent temperature and 
excitation energy for SMM and MMMC filtered predictions respectively, evaluated with 
eq.s (D and (11). 
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Fig. 9: 

Two-fragment correlation functions of the reduced velocity 1 + R{vred) for 4 < Z < 15. 
Panel a) shows MMMC calculations and panel b) SMM calculations. Solid points represent 
calculations at e* = 4 A.MeV, circles at 6 A.MeV. 
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Fig. 10: 

Experimental two-fragment correlation functions of the reduced velocity 1 + R{vred) for 
4 < Z < 15 and for e* < 3.5 A.MeV (upper panel) and for 5 < e* < 7A.MeV (lower 
panel). 
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Fig. 11: 

Partition temperature vs. e* logarithmic scatter plot. Upper panel refers to experimental 
data, lower panel to SMM predictions. The size of the squares is proportional to the yield. 
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Fig. 12: 

Experimental effective Cy as a function of the mean partition temperature. 

The calculations have been made at a density po/3 in the SMM-like hypothesis (top row) 

and at a density po/6 in the MMMC-like hypothesis (bottom row). 

The collective energy has been subtracted from the kinetic energy of the detected fragments 

in panels a), c), e) and g) and it has been included in panels b), d), f) and h). 

The level density has been taken a = f{A) in the first and second columns and A/8 in the 

third and fourth columns. 
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Fig. 13: 

Effective Cy as a function of the mean partition temperature for SMM events. 
The sohd line corresponds to e* < 25 A. MeV (flat distribution), the dashed line to 
e* < 8 A.MeV (flat distribution), the dotted line to e* < 8 A. MeV with an exponen- 
tial probability distribution. 
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Fig. 14: 

a) Experimental Campi scatter plot. The logarithm of the size of the largest fragment 
ln{Zf)ig) is plotted versus the logarithm of the second moment ln{m2) (normalized to the 
charge of the source). Three cuts are used to select the upper branch (Cut 1), the lower 
branch (Cut 3) and the central region (Cut 2). 

b) Experimental (0,e*) correlation in Cut 1 (blue), 2 (red) and 3 (green). 

c) SMM Campi scatter plot for filtered events. 

d) Partition temperature-e* correlation in the three cuts for SMM filtered events. 

The black line in panels b), d) represents the mean correlation for all the events. The size 
of the squares is proportional to the yield. 
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Fig. 15: 

a) In(Ziyig) versus ln{S2)- Fission events have been removed. The lines represent the fit of 
each branch. 

b) 5*3 versus 5*2. The hne is the fit of the correlation. 

c) 5*3 — 5*2 for the events falling in regions 1 (squares), 2 (points) and 3 (open circles) of 
the Campi scatter plot. The line is the fit of the whole correlation, shown in panel b). 

d) Charge distribution for events belonging to the Cut 2 of the Campi scatter plot. The 
line represents the power-law fit of the distribution. 
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Fig. 16: 

a) Zhig vs. e* correlation, b) ln{Z}jig) vs. /n(e*^jj — e*), c) Z^jg vs. Tcrit correlation, d) 
ln{Zbig) vs. ln{Tcrit -T). 

The solid symbols represent the values used for the fits, the lines are the fits resulting from 
equation (^. 
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